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1. Introduction 

In these proceedings we discuss continuum physical results for light hadronic quantities ob- 
tained through a combined continuum extrapolation of the RBC and UKQCD collaboration's 
32 x 64 and 24 3 x 64 domain wall fermion lattices with L s = 16 and the Iwasaki gauge action 
at f3 = 2.13 and 2.25 respectively. The lattice spacings, as determined by the combined analysis, 
are around 2.32 GeV and 1.73 GeV for the 32 3 x 64 and 24 3 x 64 lattices respectively, such that 
the lightest unitary pion masses are around 295 MeV and 330 MeV. 

The layout of these proceedings is as follows: We first discuss scaling between the two lattices, 
matching at masses within the range of our simulations. We then discuss our simultaneous con- 
tinuum and chiral extrapolation, determining the physical quark masses, lattice spacings and the 
continuum limit of the fitted observables. We make use of reweighting [[II]] in the strange sea sector 
with a corresponding interpolation in the valence sector to reach a physical unitary strange quark 
mass. We then determine the neutral kaon mixing amplitude for which we quote a preliminary 
physical value. 

2. Scaling at simulated mass points 

We compare our two lattices by matching dimensionless quantities at a simulated (unphysical) 
mass point. The method is as follows, (i) we choose a match point on either lattice (Lmatch), 
defined as a choice of simulated light quark mass m/ and heavy quark mass (ii) we interpolate 
two dimensionless quantities R\ and Rf, on the other lattice (L; nterp ) until they match those on L matc h. 
These quantities are chosen to be ratios of lattice quantities which are sensitive to the desired quark 
mass. Examples include /?/ = Mpsji/Mq and Rf, = Mps^i /Mq, where Mpsji and Mpsj,i are the 
unitary light-light and heavy-light pseudoscalar (PS) masses respectively. Prior to the interpolation, 
the numerator and denominator quantities of both /?/ and Rf, are fit linearly in nif and m/,, and the 
interpolation is performed by varying these parameters until the quantities match. 

Once the match point has been found, we can determine the ratio of the lattice spacings Z a by 
taking the ratio of a lattice quantity Q a (e.g. ciMq) between the two lattices Z a = ^Hf mteip ! . We can 
also determine the ratio of the physical quark masses 

Z — m r yS (^interp) _ 1 "M^interp) ^ _ "^ hyS (L int erp ) _ 1 % (Ljnterp ) ^ 

1 mf ys (L match ) Z u mi(L match ) ' h mf ys (L match ) Z a m h (L match ) ' 
where m = m + m res . Figure [l] shows Zj, Z/, and Z c , for various choices of match point, /?///, and Q a , 
showing good consistency between the choices. 

Picking a representative Z/ and Z/, we can determine the matching quark masses for any given 
simulated point on L matc h. We compare dimensionless ratios of other quantities over a number of 
match points in order to demonstrate the scaling of our data. Figure ^| shows two match points on 
the 24 3 x 64 lattice, demonstrating scaling violations of less than 2%. Here the Wilson parameter 
ro was obtained as per ref. [§]. 

3. Simultaneous continuum fits 

In order to make the best use of our data, we perform simultaneous fits to both lattices. We 
fit all quantities in physical units in a fixed renormalisation scheme. The renormalisation constants 
for the light and heavy quarks can be left as fit parameters or can be frozen to a predetermined 
value. For this analysis we assume the light and heavy quark renormalisation constants are equal, 
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Figure 1: Plots of Z/, Z/, (top) and Z a (bottom) determined at several match points labelled A-E on both en- 
sembles, using differing ratios R/ and /?/, and Q a (labelled in order). In the right-hand panes, for which L ma tch 
and Linterp are swapped with respect to the left, we plot the inverse quantity Z7,l , for easier comparison. 
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Figure 2: A comparison of different dimensionless quantities Q between the lattices at two match points 
(m h m h ) = (0.005,0.04) (left) and (0.01,0.04) (right) on the 24 3 x 64 lattice. Z, and Z h are calculated using 
Ri = roMps u and Rf, = r^Mpshl- 632/624 being unity defines the line of perfect scaling. 

which is justified by the results shown in figure [j]. This is also expected from renormalised field 
theory, which predicts that a single fermion mass rescaling is needed to relate theories with different 

cutoff scales. We fix Z32 = 1 and Z24 = ' hvs — r = Zi from the scaling analysis of the previous 

m, (24 3 ) 

section. This is equivalent to a non-canonical choice of renormalisation scale which can of course 
be corrected to MS at a later stage. The lattice spacings (and physical quark masses) are self- 
consistently determined by fixing a set of quantities (detailed below) in the continuum limit to their 
known physical values within an iterative procedure. 

We simultaneously perform a partially-quenched (PQ) fit to the following quantities: fps tX y> 
fps,hy, Mps^xy, Mpgfiy and Here m x and m y are PQ light valence quark masses. Following our 
standard procedure [Q] we include only light PS data for which the corresponding pion mass is 
less than ~ 420 MeV. In order to allow for a continuum extrapolation, the fit forms are obtained 
through a double expansion in ChPT and a 2 , taking terms of order m q a 2 as higher order such that 
only the LO LEC terms have a lattice spacing dependence. The PS quantities are fit to NLO SU (2) 
PQChPT fit forms and Mq is fit to a form linear in the light quark mass. For example, the NLO 
SU (2) PQChPT fit form for the light PS decay constant in the case of degenerate valence quarks is 



fps r 



f(l+c f a 2 } 



W7T g ^aT + 7 4X1 + l 5Xx ■ 



(3.1) 
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Figure 3: Results for Mps,xy (left) and fps,xy (right) on the 32 3 x 64 m/ = 0.004 ensemble obtained from 
global fits using NLO PQChPT (upper) and leading-order, linear analytic expansion (lower) fit forms. Un- 
filled symbols indicate data that is not included in the fit. The curves show the fit form at the parameters equal 
to the data set of the same colour. The red curve is the unitary mass curve at which m x = m y = mi. Note that 
the curvature apparent in the lower-left pane is a consequence of plotting 2Mp S / (m x + m y ) against only a 
single term of the denominator; m y . 

where % q = 2BZim q . The fits are performed at the physical as determined self-consistently 
through an iterative procedure with the physical kaon mass as input, using reweighting to iterate 
the sea strange mass while interpolating in the valence sector. The physical mi is obtained by 
inverting the pion fit form on the physical pion mass. The two lattice spacings are similarly fixed 
by requiring that the predicted continuum Mq matches its experimental value and that the ratio of 
lattice spacings is equal to that determined from the scaling analysis of the previous section. 



Figs. |3(a)| and [3(b)| show the NLO PQChPT fits for Mps tXy and fps. xy , not including finite 
volume (FV) corrections. Extrapolating to the continuum limit, we predict f n = 0.119(3) GeV 
and fx = 0.147(3) GeV. Comparing to the known physical values we see that these are ~ 8% 
(~ 3a) and ~ 5% (~ 2a) too low respectively. The inclusion of FV effects on f„ obtained using 
FVPQChPT |0j, provides an upwards shift by ~ 2.5%, which is insufficient to account for the 
difference. However, this disparity is of the order expected for the NNLO ChPT contributions 
obtained by squaring the NLO contribution, suggesting that a full NNLO fit is required to correctly 
reproduce the physical point. An investigation of the full NNLO fit is detailed by Mawhinney [ffQ], 
with the conclusion that these fits are not stable when applied to our data. 

Although our final fits are uncorrected, we also performed correlated fits using covariance 
matrices estimated by including increasing numbers eigenvectors Jl2| ] finding no significant change 
in our results within our limited ability to estimate the covariance matrix. 

We also detail an investigation using a form of analytic expansion about a non-zero unphysical 
pion mass, as advocated by Lellouch [ph. Using this approach we lose the ability to take the chiral 



limit and only extrapolate to the non-zero physical point. In figs. |3(c)| and |3(d)| we show the results 
for Mps.xy and fpsjy obtained using only the leading order (linear) terms of the analytic expansion. 
As an example, the LO analytic form for fps )X y is 

fps,xy = /c(l + f a a 2 ) + fiZimi + f v Zi(m x + m y ) /2 . (3.2) 
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Figure 4: Plots of our f K (left) and Bk (right) data on both lattices corrected to the continuum limit and 
the physical strange quark mass, overlayed with the NLO ChPT fit results (orange) and the LO analytic fit 
results (grey). On the former the physical value of f n is marked by a horizontal bar. 
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Figure 5: A comparison of simultaneous SU (2) PQChPT fits (upper) and LO analytic fits (lower) to the Bk 
data on the 32 3 (left) and 24 3 (right) lattices reweighted to the physical strange quark mass. 



These fit forms provide a better description of the data in our mass range than the NLO ChPT forms 
(c.f. figs. 3(a)| and |3(b)| ) even to masses significantly higher than the 420 MeV chiral cutoff. Figure 



4(a)| shows a comparison between the two fit forms (without FV corrections) on the approach to the 
continuum limit. The LO analytic fit gives predicts f n = 0.129(3) GeV and f K = 0.151(3) GeV 
which agree very well with the known physical values. 

We note that fitting an analytic expansion to higher masses and continuing this down into a 
region with known chiral non-analytic terms is not without risk, even though it worked surprisingly 
well for f % . With this in mind we continue to take the central value from NLO ChPT and use these 
analytic fits to estimate a robust error: Based on the observation that the NLO corrections to f n are 
negative in sign and assuming the convergence of ChPT we conclude that the LO analytic fits are 
likely a good method for estimating the upper bound on the systematic errors for quantities such as 
Bk for which there is no experimental value. Results at lighter masses are surely required to bring 
down the errors estimated by the spread of these two methods. 

4. The continuum limit of Bk 

The neutral kaon mixing amplitude Bk is defined as Bk = ^t^^ , where the four-quark 
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Source 


24 J Mag. 


% err 


New Mag. 


% err. 


stat 


0.010 


1.9 


0.006 


1.1 


ChPT 


0.010 


1.9 


0.013 


2.4 


FV 


0.005 


1.0 


?? 


?? 


NPR 


0.013 


2.5 


0.013 


2.4 


Scaling 


0.021 


4.0 





0.0 


Unphys. m/j 


0.005 


1.0 





0.0 


Total 


0.030 


5.7 


0.019 


3.5 



Table 1: A breakdown of the quoted error on Bk (right) compared to that of our previous determination [|4j 
(left). 

operator ffyv+AA is defined as &vv+aa = {sj^^y^d) + (sy^y^^y^d). Using the operator 
product expansion, Bk can be perturbatively related to £k, the measure of indirect CP-violation 
in the kaon sector[[I]]. The relationship contains the CKM matrix phase 8 which parameterises 
CP-violation in the Standard Model. Sk is known precisely from measurements of Kl — > 7T7T and 
Kg —>■ Tin decays, thus the measurement of Bk places constraints on the CP-violating phase. 
We fit B K on both lattices simultaneously to the SU(2) PQChPT fit form 



Bk — Bl 



l + c a a +—zr + ' 



Xy i [ Xy 
32k 2 f 2 g \A 2 



(4.1) 



where the kaon has been coupled into the theory at leading order in the relativistic expansion. Here 
Co and c\ are LECs and again we have included a term c a parameterising the &{a 2 ) corrections. 
We simplify the fit by freezing the LECs B (in % q ) and / to the values obtained from the global fit 
analysis. As discussed in the previous section, we also perform a leading order analytic expansion 
fit to the data in order to estimate the systematic error on the fit function. The fits are performed 
to renormalised lattice quantities. We determine Zb k using multiple non-exceptional momentum 
renormalisation conditions, which suppress chirality mixing. We used volume sources for good 
statistical precision. The conversion to MS has been calculated by RBC-UKQCD to one loop for 
each of the schemes, and comparison of these gives a robust estimate of higher order errors (c.f. 
refs. |f] and [g]). 

Figure |5| shows a comparison between the ChPT and analytic fits to our data at the physical 
strange quark mass. Here, the data on both lattices appears to display a slight non-monotonicity in 
the light quark mass, which on the 24 3 lattices reverses the mass-ordering of the data with respect 
to that on the 32 3 . This effect is sub-statistical and we expect it to be resolved by an increase 
in the number of measurements. Figure 4(b) shows a comparison between the unitary (m y = m/) 
extrapolations of the data in the continuum limit using ChPT and LO analytic fit forms. This figure 
encapsulates our philosophy that the extrapolation may turn down or it may not, we simply do 
not know but we ensure that our errors cover this difference. Here the difference between the two 
physical points is much smaller than that for f n , around 2.4%, but is again of a similar magnitude 
to the expected NNLO ChPT contributions. 

Using the lattice spacings and quark masses obtained from our global fit procedure we quote a 
preliminary physical value of B^ s (2GeV) = 0.537(19), not including finite-volume effects. These 
results will be published in our forthcoming paper |[T3|]. Comparing this to our previous value of 
B^ s (2GeV) = 0.524(30)0] we see a consistent result with a substantially improved total error. 
Table |I] contains a breakdown of the total error comparing to our previous result. Note that with 
the inclusion of reweighting in the heavy quark mass and by including iff (a 2 ) corrections in our fit 
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forms, systematic errors associated with the unphysical mj, and with scaling violations have been 
removed and replaced with contributions to the statistical error. 
5. Conclusions and Acknowledgements 

By matching the data on our 24 3 x 64 and 32 3 x 64 lattices at an unphysical match point, we 
have shown that our data is scaling to within 2% for all quantities we tested. With confidence in our 
data we then described a global fit procedure whereby the quantities Mpsji, fps.n, Mps,hh fps,hl an d 
Mq are all fit simultaneously on both lattices to NLO ChPT fit forms including a parameterisation of 
the leading (a 2 ) effects. This global fit is iterated within a procedure that determines the physical 
quark masses and lattice spacings. We make use of reweighting in the strange sea quark mass 
in order to reach the physical point. We show that NLO ChPT fit forms do not reproduce the 
known continuum value of f K , giving a value around 5% (~ 2a sta t) too low after finite volume 
effects are included. We demonstrate that an analytic expansion (as advocated by Lellouch [§]) 
reproduces the physical f % and provides a better fit to the data at mass scales above the chosen chiral 
cutoff. Noting the risk of extrapolating these analytic forms down into a region with known chiral 
non-analyticities, we choose to continue to take the NLO ChPT central value for the subsequent 
analysis of the neutral kaon mixing parameter Bk, taking the LO analytic fits to this quantity as an 
upper bound on the systematic errors. We emphasise that results at lighter masses are necessary 
to reduce these errors. Finally we give a preliminary continuum result for Bk as B^ s (2GeV) = 
0.537(19), where the 3.5% quoted error includes all systematics bar the finite-volume corrections 
which have not yet been calculated. We use the renormalisation parameter Zb k calculated using 
multiple non-exceptional momentum schemes with volume sources for good statistical precision. 
The MS conversion factors have all been calculated to one loop giving a robust systematic error 
estimate. We compare our Bk value to our previous published result 

We thank all members of the UKQCD and RBC collaboration. Computations were performed 
on the QCDOC machines at the University of Edinburgh and Columbia University, the US DOE 
and RBRC facilities at the Brookhaven National Laboratory and at the Argonne Leadership Class 
Facility. The resources made available at the Argonne Leadership Class Facility were essential to 
the generation of our 32 3 x 64 ensembles. The author was supported by the UK STFC. 
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